Critical States in a Dissipative Sandpile Model 
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A directed dissipative sandpile model is studied in the two- dimension. Numerical results indicate 
that the long time steady states of this model are critical when grains are dropped only at the top 
or, everywhere. The critical behaviour is mean-field like. We discuss the role of infinite avalanches 
of dissipative models in periodic systems in determining the critical behaviour of same models in 
open systems. 
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Spontaneous emergence of long range spatio-temporal 
correlations in driven dynamical systems without fine 
tuning of any control parameter, is the concept of Self- 
Organized Criticality (SOC) Since its introduction 
m 1987 0, the precise conditions which are necessary and 
sufficient for SOC, are subjected to intense scrutiny. The 
question which attracted much attention is, can one have 
criticality if there is a non-zero rate of bulk dissipation? 
While some works at the early stages Q] suggested that 
indeed, the conservation of the transported quantity in 
the dynamical rules is a necessity, the later works claimed 
a negative answer. 

In this paper, we study a directed dissipative sandpile 
model and our numerical results indicate that it is crit- 
ical. We argue, that a dissipative model may be critical 
provided the dissipation is not too strong and conjecture 
a criterion to determine the critical behaviour. 

In the sandpile model of SOC, sand grains are locally 
injected and transported on an arbitrary lattice. Too 
many grains cannot be accommodated at any site. A 
site relaxes if the number of grains exceeds certain cut- 
off and transfers the grains equally to the neighbouring 
sites. This transfer process is conservative, where no 
grain is lost or created. At the critical state, cascades 
of relaxations follow due to single injection of grains, 
which are called avalanches. Grains, however, dissipate 
out of the system through the boundary, otherwise no 
steady state is possible. This is called the Abelian Sand- 
pile Model (ASM) |],[| . A globally driven conservative 
earthquake model is also similarly defined where energy 
is fed uniformly at all sites and transported ||. This 
model reproduces power laws of energy release similar to 
the Gutenberg- Richter law || . 

There are some studies on the dissipative models also. 
MKK studied a sandpile model where a grain can dis- 
sipate during a relaxing event, in a probabilistic man- 
ner. Numerical findings show that the system reaches 
a sub-critical state with the characteristic sizes of the 
avalanches depending inversely on the probability of 
dissipation |l(J. On the other hand, the dissipative 
ASM showed criticality with mean-field like critical be- 



haviour jiT| . One dimensional version of this model also 
showed critical behaviour even with finite driving rate 
Q . The OFC model JTJ] studied the dissipative earth- 
quake model, where dissipation is controlled by a pa- 
rameter a. It is claimed that the OFC model is critical 
for a c < a < a , the conservative value of a being a , 
with the critical behaviour depending on a Jl3|-[l6| . The 
stochastic version of OFC model, however, is shown to 
loose criticality for any a < a pj\j. 
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FIG. 1. Configurations of DDSM on a lattice of size 
L = 512. (a) Avalanches created by dropping grains ran- 
domly only at the top row, orderly place grains along 'V 
shaped lines, (b) Grains are dropped at all sites. Different 
densities in different regions indicate the age of the big aval- 
nches passed through that region. 



In a conservative model of SOC the grains move a dis- 
tance of the order of the system size L when started 
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from the inner most region. This makes the average 
avalanche size grow as a power of L, so that an infi- 
nite system has a power law distribution of the avalanche 
sizes. In contrast, in a dissipative model, the grains dis- 
sipate at any distance within the system. If all grains 
do dissipate within certain cut-off distance, the average 
avalanche size would not have any dependence on L in 
the large limit. Therefore, for a dissipative model to be 
critical, only a fraction f(L) of grains should dissipate 
from the bulk and the rest through the boundary, such 
that, /(oo) = _Lti^oo/(_L) < 1. We present examples of 
both cases in the following. 



variable length since the fluctuations take place only at 
the two ends. If the TL has a length £ at time t, it would 
have the length at least £ — 1 at time (t + 1), since all 
inner (£—1) sites will get two grains each and will topple 
again. If either or both neighbours of the end sites are 
occupied, the TL collects the grains into it, and grows 
in length to £ or, i + 1. However, if these neighbours 
are vacant, TL fills them and shrinks in length. There- 
fore, the two ends of the TL move in principle like two 
annihilating random walks starting from the same point. 
The avalanche terminates when they meet and annihilate 
each other Q . 
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FIG. 2. Plots of P(t)t Tt ~ x (circles), p(y)y a (squares) and 
(j(y)y 13 (diamonds) for system sizes L = 2048 (filled sym- 
bols) and 13000 (opaque symbols). Horizontal portions of the 
curves correspond to t% = 2.027, a = 1.012 and j3 = 1.184. 



On an oriented square lattice with extension L, sites 
are either vacant or occupied by single grains in the stable 
state. System is periodic along the x direction and the y 
coordinate increases downward. Grains are dropped ran- 
domly. A toppling occurs only when the number of grains 
hi > 1, the site i is then vacated: hi — > 0. The system 
has a preferrence along the y direction and the down-left 
and the down-right neighbours at the next row gets one 
grain each: hj —>hj-\-l. In a toppling with height 2, 
grain number is conserved, where as, in a toppling with 
height 3, one grain dissipates from the system. Unlike 
the Directed Abelian Sandpile model (DASM) Q, our 
model is non-Abelian since sites are vacated in a top- 
pling and we call it as the 'Directed Dissipative Sandpile 
Model' (DDSM). 

In a parallel dynamics all toppling sites reside in a 
single row on a contiguous toppling line (TL). It has a 
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FIG. 3. Saturation of the density p(y) and the fraction of 
bulk dissipation a(y) in a system of size L — 13000. 



Grains are randomly dropped in two ways: In case A, 
they are dropped only on the top row at y = 1, and in 
case B, they are dropped everywhere. Grains dissipate 
through the boundary at y = L. First, we consider case A 
and a stable configuration is shown in Fig. 1(a). Grains 
are marked by black dots, where as vacant sites are made 
blank. Lines of grains in the shape of 'V are mostly ob- 
served. This is because, due to the bulk dissipation, the 
density is so low that the TL moves almost in a deter- 
ministic manner. A 'V is formed by the movement of 
a TL through a vacant region. In this case the TL uni- 
formly shrinks, leaving behind a trail of two converging 
lines of occupied sites at the two ends. However, a TL 
may also propagate in a 'A' between two Vs. In that case 
it uniformly grows in length, deletes two sides of two Vs 
up to their lowest points and then starts shrinking, pro- 
ducing two converging lines which finally make a bigger 
V. In this way bigger V shapes are generated at the ex- 
pense of smaller Vs, which finally reaches the boundary 
at the bottom and dissipates. Such almost determinis- 
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tic dynamics makes avalanches of rectangular shape in 
general, but mostly they are squares. 

An avalanche deletes all occupied sites through which 
it passes. No dissipation occurs on the first two rows 
where the average density is 1/2 as in DASM fjlf]. It 
then decreases with y as a power law: p(y) ~ y~ a . A 
system of size L = 13000 is simulated by dropping 2 x 10 9 
grains. We plot p(y)y 1012 with y on a double logarithmic 
scale, the curve is horizontal for the large y values, giving 
a = 1.012 ±0.030 (Fig. 2). 




FIG. 4. The avalanche size distribution for the case B of 
DDSM. The dot-dashed curve is for the high density region. 
The four curves with solid lines are for avalanches of the whole 
system for the system sizes L = 256, 1024, 4096 and 13000 
(from top to bottom). 



Avalanche size s is the number of sites toppled in an 
avalanche. Simulation results indicate that the cumu- 
lative probability distribution of s follows a power law: 
P{s) ~ s 1_T ' with t s = 1.52 ±0.03. The life-time t of an 
avalanche is its vertical extension along the preferred di- 
rection and also follows similar power law: P{t) ~ t x ~ Tt 
with r t = 2.027 ±0.030 (Fig. 2). The average avalanche 
size < s(t) > varies with life-time t as < s(t) >^ t lst 
with j st = 2.01 ± 0.03. Since s, t are two measures of 
the same random avalanche cluster, they are necessar- 
ily dependent, and are related by the scaling relation: 
7st = (r t - l)/(r s - 1). 

We explain these results in the following way. It is 
reasonable to assume that most of the avalanches are of 
rectangular shapes, which implies that j s t = 2. Now, if 
the TL has a width w(t') at the intermediate time t', then 
2w(t') grains cross that row y = t' . The dissipation flux 
per grain can be divided into 'bulk-flux' and 'boundary- 
flux'. All grains crossed by the TL, except at its two ends, 



dissipate. Therefore, the density and the system size L 
controll the share between the bulk and the boundary 
fluxes. The constant average boundary-flux through the 
row at y is < w{y) > y 1 ~ n which gives < w(y) >^ 
y Tt ~ x . But, since average avalanche size of life-time t is 
< s(t) >= ^w(t')dt' = t Tt , we get j st = r t = 2 and 
t s =3/2. We numerically check the relation < w(y) >= 
ky and find a nice straight line with slope k — 0.312 ± 
0.001 and the correlation coefficient 0.999. 
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FIG. 5. (a) ASM on a periodic 2x2 lattice, which leads to 
a periodic infinite avalanche (b) Dissipative ASM on a similar 
periodic lattice also leads to the periodic infinite avalanche. 



Since the density of the system decreases with increas- 
ing y, we expect that the average dissipation also should 
decrease with increasing y. The fraction cr(y) of total 
number of grains dissipated in the y-th row varies as: 
v{y) ~ V' 13 , where (3 = 1.184 ±0.030 (Fig. 2). Therefore 
the bulk-flux f(L) should vary as: f(L) — f(oo) — CL~ X 
with x = (3 — 1. The exponents x is estimated inde- 
pendently by plotting f(L) vs. L~ x for different x val- 
ues and the best value obtained is x = 0.17 ± 0.03 and 
/(oo) = 0.634 ±0.010. 

Now we consider the case B (Fig. 1(b)). The local 
density fluctuates widely since an avalanche sweeps the 
region it passes and the local density in this region re- 
starts growing from the scratch. Since the grains are 
dropped everywhere uniformly and the bulk dissipation 
depends on the density, we expect that there should be a 
saturation region where the average density is constant. 
The density p(y) decreases from 1/2 and then saturates 
to a constant value 0.1543±0.0010 around y c w 100. The 
rate of dissipation cr(y), initially increases but finally sat- 
urates to the uniform dissipation limit: a(y) = C/L, with 
C = 1.01 (Fig. 3). The bulk-flux f(L) asymptotically 
reaches to /(oo) = l as 1/L. 

It turned out that the system has two regions. The 
high density region extends from the top to y c and the 
saturation region from y c to L. We separately collect the 
distribution data for the avalanches originated in these 
two regions. For the high density region, the t s « 1.5 and 
T t w 2.0 are obtained as in the case A and < s(L) >~ 
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L and < t(L) >~ logL are observed. Linearity in < 
w(t) >= k\t is still obeyed with fci=0.1 and j st ~ 2 
is obtained again. However, for the saturation region, 
plots of the distribution data showed two regions: an 
initial high slope 2.5 for the small s values, followed 
by a slope t] « 1.5 for the large s values. We explain 
the large value of is due to those small avalanches, 
which grow on an empty region swept out by a previous 
large avalanche. When this region reaches the steady 
state, the avalanches get the usual exponent t 1 s = 1.5 
for large s values. We see that both < s(L) > and < 
t(L) > have constant values independent of L. The total 
distribution, has the behaviour of the saturated regions, 
since the avalanches generated in this region have larger 
weights. 

We now look into the effect of the boundary on dissi- 
pative models in more detail. In a conservative sandpilc 
model with periodic boundary condition, the total mass 
of the system grows up indefinitely. Very soon, an 'Infi- 
nite Avalanche' starts which never terminates. For ASM 
on a periodic square lattice, the same height configura- 
tion repeats at certain interval, toppling all sites exactly 
once. The period is of the order of L and is dependent 
on the initial configuration. We show such a 2 x 2 system 
in Fig. 5(a). Next we consider the dissipative ASM ]Tl| ] 
on the same lattice. After some initial dissipation, this 
model also creates a periodic infinite avalanche which is 
dissipationless (Fig. 5(b)). We now test our DDSM on 
periodic system, by making the y direction also periodic. 
We observe again dissipationless infinite avalanches in 
both cases A and B. A TL in the form of a ring moves 
indefinitely with uniform speed on the empty torus. 

An infinite avalanche has to be dissipationless after 
some time, otherwise it will make the whole system 
empty. If a dissipative model in a periodic system has 
no infinite avalanche, it indicates that even the largest 
avalanche is not as big as the system. Therefore for the 
same dissipative model in the open system, the bound- 
ary should have no effect on the avalanche sizes, leading 
to sub-critical states. We conjecture that: A dissipative 
model will not show self-organized criticality if the same 
model on a periodic system has no infinite avalanches. 

To verify this conjecture, we check some examples. The 
probabilistic dissipation model pOf , the stochastic OFC 
[17 1 model and the dissipative two states sandpile model 
[19 1 all are non-critical on open boundary systems and 
do not have infinite avalanches on the periodic systems. 
However, the random creation-dissipation model in [ jTo) , 
the dissipative ASM [jll| , the cases of DDSM as described 
in this paper, lead to SOC states with open boundary and 
also have periodic avalanches on the periodic systems. 
Finally, we check that the deterministic OFC model |l3| ] 
also does not produce any infinite avalanche on the pe- 
riodic system for any a < 1/4. Therefore, according to 
our conjecture, deterministic OFC model is not critical, 
which is against the general belief. In a recent preprint, 



it has been claimed that the deterministic OFC model is 
critical in the conservative regime only ppfl . 
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